Nucleon spin structure and mass of quarks 



Petr Zavada 



Institute of Physics, Academy of Sciences of Czech Republic 
Na Slovance 2, CZ-180 40 Prague 8 
E-mail: zavada@fzu.cz 

May 31, 1999 

Abstract 

The alternative to the standard formulation of the quark-parton model is pro- 
posed. Our relativistically covariant approach is based on the solution of the master 
equations relating the structure and distribution functions, which consistently takes 
into account the intrinsic quark motion - in contradistinction to the standard infinite 
momentum approach, in which this motion is latently suppressed. The model well 
reproduces the experimental data on the both polarized and unpolarized structure 
functions, assuming that only the valence quarks term contributes to the nucleon 
spin. It is shown, the combined analysis of the polarized and unpolarized data can 
give an information about the effective masses and intrinsic motion of the quarks 
inside the nucleon. Simultaneously, it is shown that the rate of the nucleon energy 
carried by the quarks can be less, than estimated from the standard approach. As 
an addition, a prediction for the proton spin function c/2 is given. 



1 Introduction 

Measuring of the structure functions is an unique tool for the study of the nucleon internal 
structure - together with the quark-parton model (QPM) giving the relations between 
the structure functions and the parton distributions, which represent a final, detailed 
picture of the nucleon. In this sense, these relations, obtained under definite assumptions, 
are extraordinary important, since the distribution functions themselves are not directly 
measurable. At the same time, the standard, simple formulas relating the structure and 
distribution functions are ordinarily considered so self-evident, that in some statements, 
the both are identified. 

The experiments dedicated to the Deep Inelastic Scattering (DIS), are oriented to 
the measuring either unpolarized or polarized structure functions. The results on the 
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unpolarized functions are compatible with our expectations based on the QPM and QCD, 
but the situation for the polarized functions is much more complicated. Until now, it is not 
well understood, why the integral of the proton spin structure function g 1 is substantially 
less, than expected from very natural but possibly equally naive assumption, that the 
nucleon spin is created by the valence quarks. Presently, there is a strong tendency to 
explain the missing part of the nucleon spin as a result of the considerable contributions 
of the sea quarks (particularly strange quarks) and the gluons. Nevertheless, a consistent 
explanation of the underlying mechanism is still missing. During the last years, the 
hundreds of papers have been devoted to the nucleon spin structure, for the present 
status see e.g. [[!]], 0, the comprehensive overviews and citations therein. 

In the present paper we continue our discussion started in ||, ||, where we have shown, 
that the standard formulation of the QPM, conceptually firmly connected with the infinite 
momentum frame (IMF), oversimplifies the parton kinematics. In we demonstrated 
that the effect of oversimplified kinematics in IMF can have an impact particularly on 
the spin structure function g 1: or more exactly, it can substantially modify the relation 
between the distribution and spin structure functions. We have suggested this effect 
can be a source of the discrepancy between the experimentally measured and naturally 
expected magnitude of the spin function g±. The primary aim of this paper is to precise 
just this point. 

The paper is organized as follows. In the following section we add some comments 
to our master equations connecting the structure and distribution functions and point 
out the distinction between them and these used in the standard approach. In the Sec. 
H]. we propose the model, in which the internal motion of quarks is consistently taken 
into account. Further, in contradistinction to the standard treatment based on the QCD 
evolution of the distribution functions in dependence on Q 2 , our model rather aims at 
describing the part of distribution and structure functions, which are not calculable in 
terms of the pertubative QCD. In the Sec. (|. the results of the model on the polarized 
and unpolarized proton structure functions are compared with the experimental data 
and some free parameters are fixed. Some important additional comments on the model 
and obtained results are done in the Sec. ^. The last section is devoted to the overall 
summary and concluding remarks. Since this paper should be read together with 0, ||, 
for convenience we refer to the equations or figures in these papers simply with prefixes 



For completeness, in the Appendix we give a more rigorous explanation, that the 
approximation of the term 



entering the argument of the 5— functions in some integrals in this and two the previous 
papers is for our purpose fully acceptable. 
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2 Master equations 



In the preceding discussion @,@ we have shown [see Eqs. (P3.41), (Q2.1)], that if 
one assumes momenta distributions (Q2.2), (Q2.3) of quasi-free quarks having mass m 
are spherically symmetric in the nucleon rest frame, then the corresponding structure 
functions W±, W 2 , Gi, G 2 obey the master equation 



W 2 \ 



s a MG 1 + (Pqs a -sqP cr )^ 



+ A(P a qp + P p q a ) + Bq a qp 
cPp 

G{po)(2p a p p - g a/ 3Pq)S((p + q) 2 - m 2 )- 



Po 

+ie a px«q X I H{pv)mw a 5{{p + q) 2 -m 2 )^. (2.1) 
J Po 

Simultaneously, we have shown, how the explicit solution of this equation looks like. In 
this solution, the mass xq = m/M formally appears as a free parameter. Let us make a 
remark, in the same way a similar equation can be obtained and solved also for the set 
of the neutrino structure functions, nevertheless in this paper we shall deal only with the 
electromagnetic ones. 

Let us note, despite of the fact that Eq. (|2.1|) is assembled for quark momenta dis- 
tributions G, H in the nucleon rest frame, the equation is relativistically covariant. Its 



manifestly covariant form follows immediately from (2.1) after the substitution 



G{Po Mb ) = G , H(p 0;lab ) = H (^0 . (2.2) 



For moving nucleon we have 



P = (P ,P 1 ,0,0), 7 = 5 (2-3) 



and 

pP _ PqPq - pxPi 
M ~ M 



P ' M 



lipo ~ PPl) = Po,lab (2.4) 



which means, that the phase space of the subset of quarks with po,Za& fixed, represented 
by the sphere 

Pi+pl+pl=pl,iab- m2 (2-5) 

in the nucleon rest frame, is in a boosted system correctly represented by the ellipsoid 
with the shape defined by the Lorentz transform ( p.4| ). 
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Further, the Eq. Q2.1D involve also the quark polarization vector w^, which as follows 
from 0, plays the crucial role for the evaluation of the spin structure functions gi,g 2 - So 
here we consider desirable to derive its form in a more rigorous way, than we suggested for 
Eq. (Q2.12). Generally, quark polarization vector should be constructed from the proton 
momentum P, proton polarization vector s and the quark momentum p: 

w„ = AP^ + B S p + C PfM , (2.6) 

where A,B,C are invariant functions of P,s,p. Then contracting of Eq. (Q2.5) with 
Pa,qa and s a and using relations (Q2.4) give the equations 

-sqMG 2 = ^f H (^) [AM- + CpP) s(^--x) ^, (2.7) 
H 2Mv J \M J v y ' \Mv ) po K J 



pq \ d 3 p 
— x 



x ^, (2.. 
J Po 



-MGi - vG 2 = ^- I H ( ?£) (-B + Cps) s(£L-x) d X (2.9) 



Mv ) po 



Elimination of G\, G 2 gives 



H (PP) CP ( vP/M — q _\ s m_ x) fp =Q (2 1Q) 

and since P, q, s are independent, C must be zero. The remaining invariants A, B follow 
from Eq. (Q2.4), which implies 

A 2 M 2 -B 2 = -1, ApP + Bps = (2.11) 

and solution of these equations reads 

A- T _ PS B = ± PP (2 12) 

V(P P ) 2 - {ps) 2 M 2 ' V(pP) 2 - {ps) 2 M 2 ' 

So the quark polarization vector has the form 

w, = ± <fP)s>-<Ps)P* (213) 
M ^(pP) 2 - (ps) 2 M 2 

Contributions of both possible solutions [sign +(— ) means that quark spin is parallel 
(antiparallel) to the proton spin in its rest frame] are in our calculation taken into account 
by the difference in Eq. (Q2.3). Apparently, for the proton rest frame and polarization 
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s = (0, 1,0,0) the last equation is identical to Eq. (Q2.12). Further, obtained invariants 
A, B, C now give spin structure functions from Eqs. fl2.7Q, (|2.8|) in covariant form 



= m f f pP\ (pP)(sq) - (ps)(Pq) ^ ( pq \ d 3 p 

1 2M(Pq)(sq)J \M J y/(pP) 2 - (ps) 2 M 2 \Pq J Po ' K ' ' 



G 2 = ,™*f . /"gf^ . PS S (P*- x )*l. (2.15) 

2(Pq)(sq)J \MJ ^(pPy-(p S yM* \Pq Jpo K J 

Apparently, according to these relations the structure functions can depend also on mutual 
orientation of s and q. Of course, this dependence is more complicated, apart the factor 
sq ahead of the integrals, integration involves the terms ps and pq. This question is being 
studied and will be discussed in a separate paper. In the present paper our considerations 
will be based on the results obtained in ||, which follow from Eqs. ( |2.14| ), ( |2.15| ) applied 
in the proton rest frame for the longitudinal polarization s = (0, 1,0,0). Obviously, for 
this case the last two equations are equivalent to Eqs. (Q2.13), (Q2.14). 

The scheme based on the Eqs. Q2.1J ) and (|2.13|) with all their implications suggested 
in ||,|| can be a priori valid for quasi- free quarks (on mass shell) filling up the nucleon 
volume. In this sense the scheme represents a covariant formulation of the naive QPM. 
We have shown that Eq. fl2.1|) in which the quark internal motion is consistently taken 
into account imply relations between the structure and distribution functions different 
from those obtained in the standard procedure relying on the IMF, which is based on the 
approximation p^ = x^P (parton momenta components are x fractions of the proton mo- 
menta components). In fact this approximation in the covariant formulation is equivalent 
to the assumption, that the partons are static with respect to the nucleon, therefore there 
are suppressed not only the transversal momenta, but also longitudinal ones. Of course, 
this consequence is somewhat obscured just in the IMF, where all the relative motion is 
frozen, since all the processes run infinitely slowly - including the passing of the probing 
lepton through the nucleon. Let us remark, the standard relations (e.g. F 2 = x e fqi) 
obtained in the naive QPM with static quarks are currently applied even in the standard 
approach based on QCD improved QPM, which is not a consistent procedure, since it 
means that correct dynamics is combined with incorrect kinematics. 

In this way we have shown, that the relations between the structure and distribution 
functions can be, at least on the level of the naive QPM, strongly modified (particularly 
for the polarized case) by the parton internal motion. This result can be instructive by 
itself. The impact of the quark intrinsic motion on the function gi(x) has been discussed 
also in some other approaches |7[|-||11| and necessity of the covariant formulation for the 



spin structure functions has been pointed out in |12 . 

However, the real nucleon is much more complex object, than just a bag of quasi-free 
fermions. But we shall try to assume the following. The relations obtained within the 
scheme suggested above can be used as a good approximation even for the interacting 
quarks, but provided that the term 'mass of quasi-free parton' is substituted by the term 
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'parton effective mass'. By this term we mean the mass, which a free parton would have to 
have to interact with the probing photon equally as the real, bounded one. Intuitively, this 
mass should correlate to Q 2 : a lower Q 2 allows more time and space for the struck parton 
to interact with some others, as a result the energy is transferred to a greater system than 
the parton itself. And on contrary, the higher Q 2 should mediate interaction with more 
"isolated" parton. Moreover, we accept that the value of the effective mass even for a 
fixed Q 2 can fluctuate - e.g. in a dependence on the actual QCD process accompanying 
the photon momentum transfer. This means, that the terms in the relations involving the 
mass of quasi-free parton x — ra/M will be substituted by their convolution with some 
'mass distribution' p : 



/(z„) -> J p(x ,Q 2 )f(x )dx . (2.16) 

In the following we shall propose a simple, but sufficiently general model for the 
unknown distributions p, G, H, in which all the dynamics of the system is absorbed. 
Then, these distributions will be used as an input for the calculating of the corresponding 
structure functions. 



3 Model 

Construction of the model is based on the following assumptions and considerations: 

1) Parton distribution P(e)de representing the number of quarks in the energy interval 
(e, e + de) can be formally expressed : 

P(e) = y E,r j jp j {e), J> = 1, (3.1) 
j j 

where rj is a probability that the nucleon is in the state with j charged partons (quarks + 
antiquarks) of various flavors, and pj is the corresponding average one-parton distribution, 
which fulfills 

j p 3 (e)de = I. (3.2) 

2) Nucleon consists of the three quarks and partons (gluons + qq pairs) mediating 
the interaction between them, as sketched in the Fig. [l]a, where the individual pictures 
represent some terms in the sum (|3.1|) . The flavors and spins of all the quarks in each the 
picture are mutually cancelled, up to the three quarks giving additively the corresponding 
nucleon quantum numbers. These three quarks are in the figure marked by black and in 
our approach are identified with the valence quarks. The reason, that such identification 
is quite sensible, is the following. Apparently, the sum ( |3.1| ) can be split into quark and 
antiquark parts P q (e) , Pq-(e) , then our valence term reads 

P val (e)=P q (e)-P q (e), (3.3) 
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which in the x— representation exactly corresponds to the current definition of the valence 
quarks. Correspondingly, the unmarked quarks are identified with the sea quarks. But 
both the kinds of quarks have the same energy distributions Pj(e) entering the Eq. (|3.1|) , 
in this sense they are completely equivalent. On the other hand, it is obvious, that for 
the valence quarks, in Eq. (|3.1| ) only "black" quarks from the figure contribute, therefore 
if pj is assumed in the first approximation independent on the flavor, then 



Pval(e) = 3^2r jPj (e). 

i=3 

3) The quarks carry only part of the nucleon energy (mass), 



P{t)ede = c q M, 



(3.4) 



(3.5) 



where the factor c q equals roughly one half, the rest is carried by the gluons. In the first 
approximation we shall assume this factor is valid also for any term in the sum (|3 . 1|) . 



(3.6) 



which in other words means the ratio of the total energies of quarks and gluons, together 
constituting the nucleon mass, is the same for all possible states sketched in Fig. 

4) We assume all the quarks in the nucleon state j have approximately the same 
effective mass rrij and denote Xq = rrij/M. One can expect, for higher j the parameter xq 
will drop and so the sum (|3~4T) can be substituted by the integral 



P V ai(t) = 3 J pv(x )p(e,x )dx , 



pv(x )dx = 1. 



Obviously, Eq. ( |3.1| ) can be with the use of Eq. ( |3.6| ) rewritten in a similar way: 



P(e) 



p(x )p(e,x )dx , 



where 



P{x ) = pv{x )- — - 
e(x ) 



e(x ) = / p(e,x )ede. 



(3.7) 



(3. 



(3.9) 



The physical meaning of the distributions pv,P is the following. The distribution p(xo) 
represents a probability, that the effective mass of the quark, which the probing lepton 
interacts with, is Xq or alternatively, p(xo)dxo is the number of quarks in the interval 
(xo, Xq + dxo), which the lepton has chance to interact with. On the other hand, the 
(normalized) distribution pvi x o) can be interpreted as a probability, that the exchang- 
ing photon "distinguishes" the quarks with the effective mass xq - as expressed by the 
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pictures with different granularity in Fig. |T]a. In this sense, each picture in the Fig. 

can be labeled by some xq, equally as the corresponding term p(e, Xq) in the integral 
( p.7|) . Obviously, at the same time the py(x ) represents also the distribution of effective 
masses corresponding to the valence quark term. Intuitively, the probability of different 
contributions in Fig. [I]a should depend also on Q 2 (higher Q 2 = 'better resolution'), so we 
expect 



In the next, we shall identify these distributions with that suggested on the end of the 
previous section. 

5) The relations (P3.52), (P3.56) and (P3.20) give the recipe how to obtain the struc- 
ture function F 2 from a given energy distribution of the partons with some fixed value Xq 
and charge e q : 

F 2 (x, x Q ) = e 2 q (p(x,x ), 



<p{x,x Q ) = -x 2 



1 2K(0e 

e + x 2 



M 

Y 



(3.11) 



where the sign +(— ) in the second relation refers to the region £ < xq (£ > xq). For 
the application of this procedure to Eqs. (|3.7|), ( p.8| ) one has to weight the contributions 
integrated over x$ by the corresponding (mean) charge squared. Apparently, the charge 
weight of the valence quarks is the constant 



W va l = ~ 



1 

3 



(3.12) 



for the proton and similarly for the neutron, w va i = 2/9. For the sea we assume in the 
first approximation the "equilibrated mixture" of the quarks u : d : s = 1 : 1 : 1, so 



2\ 2 /r 2 



2 
9' 



(3.13) 



Then for the nucleon with j quarks we get 

3w va i + (j - 3)w sei 



W ; 



J 



W sea + 



3(w val - W sea ) 



or in terms of Xq 



W(X ) = W sea + 3(w va l - W sea ) 



e(x ) 
c q M' 



(3.14) 



(3.15) 
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Therefore, the energy distributions (|3.7|), (|3.8| ) generate the corresponding structure func- 
tions: 

F 2v cd(x,Q 2 ) = 3w val / fi V (x ,Q 2 )ip{x,x )dx , (3.16) 



F 2 (x, Q 2 ) = j (j^^w sea + 3(w va i - w sea U /iy(x , Q 2 )^(x, x )dx . (3.17) 

6) Now, let us pay attention to the spin structure functions. According to the concept 
suggested in item 2), only valence quarks contribute to the nucleon spin. First, we shall 
consider the spin functions generated by the valence quarks with some fixed effective mass 
x , then we shall easily proceed to the case with the distribution p v (x ). 

In H , for sake of simplicity, we have assumed all the three valence quarks contribute to 
the proton spin equally [Eq. (Q2.7)]. On the other hand it is obvious the quark symmetry 
group can impose an extra constraint on the contributions of different quark flavors as it 
follows e.g. from the philosophy of the well known Bjorken |13| and Ellis- Jaffe |14| sum 
rules based on the symmetries U(6) and SU(3). In our consideration we shall not strictly 
assume any particular group of symmetry, but the different spin contributions of u- and 
d-quarks will be expressed by a free parameter a, < a < 1, having in the notation of 
Eq. (Q2.7), e.g. for the proton, the following sense 

Ah u (p ) = 2ah(p ), Ah d ( Po ) = (1 - 2a)h(p ), (3.18) 

where h is the valence distribution 

«(po) = d(p ) = %o), J h{p )d 3 p = 1, (3.19) 

which is not, due to different normalization, identical with the distribution p(e), but the 
both are simply related 



pie) = 4We 2 -m 2 h(e), (3.20) 

in the same way, as the distributions P, G in Eq. (P3.14). 
In the case of proton, there are the particular cases: 

a) a = corresponds to the mutual spin orientation of the three valence quarks (s u , s u , sj) = 
(-1,+1,+1). 

b) a = 1/3 corresponds to the oversimplified scenario studied in 0, assuming the equal 
contribution of all the three quarks; (s u , s u , Sd) = (+1/3, +1/3, +1/3). 

c) a = 2/3 corresponds to the non-relativistic SU(6) approach. From the wave function 



1/1 I 



p, t) = ^= \-y= \(ud + du)u - 2uud) ® — |(t| + |t) T -2 TT!) 
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+ _L \ {ud _ du)u) g, _L | (U _ iT) (3.21) 

one can easily show the mean value of the spin carried by the d(u)— quarks is —1/3(4/3), 

i.e. (s u , s u , Sd) = (+2/3, +2/3, —1/3), which agrees with a = 2/3 in Eq. (|3.18| ). 

d) a = 1 corresponds to the mutual orientation of the three quarks (s u , s u , Sd) = (+1, +1, — 1) 

So, the proton spin function H(p ) entering the master equation (|2.1]) and expressed 
in terms of the functions Q3.18 ) reads 



H p ( Po ) = ^au(p ) + i(l - 2a)d(p ) = ~(1 + 6a)% )- (3.22) 
Assuming the neutron is isospin symmetric, its corresponding spin function will be 

H n ( Po ) = ^(1 - 2a) M (p ) + ^ad(po) = ^(4 - 6a)% ), (3.23) 
therefore the corresponding equations for the nucleon spin structure functions read 
g?(x, x ) = vtpijjjfa x ), gj(x, x ) = w^Jj^x, x ), j = 1, 2, 

<nn = ^(1 + 6a), w n spm = i(4 - 6a), (3.24) 
where, in accordance with Eqs. (Q2.17), (Q2.18) 



Mx , n) = %[tW.fe±E S (*LgL- x \f Pl (3.25) 



2 i Po V Po - Pi V M 



m f h(p ) Pl - fpo+Pl \ ,:; 



The function (x, xo) can be with the use of Eqs. (Q2.24),(Q2.27) expanded 

M*,*o) = yE(7)(-l) 3 ^3/ 2 W (|) 5+1/2 . (3.27) 
Since Eq. (P3.56) implies 

pi / of \ i+ 3 / 2 

^•-3/ 2 (x) = - jf V '(0 (7^2 J ^ (3-28) 



one can easily show the sum in Eq. ( 3.27 ) gives 



M*,*o) = -xo / J e + ll_ x ^- (3-29) 
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Similar manipulation with the function ip 2 gives the result 



1p 2 {x,X ) = -y 



1 vzm e+xi-2xj 



(3.30) 



Obviously, for the case with the distribution /zy, the corresponding spin structure func- 
tions read 



(3.31) 



gj(x,Q ) = w spin / [j, v (x ,Q )ifj j (x,x )dx , j = 1,2. 



Let us note, the structure functions F 2 , F 2va i, g±, g% are not independent, all of them are 
in the corresponding way generated by the distributions fly and Vo (or, equivalently by 

p)- 

7) Now, to make the construction suggested above applicable for some quantitative 
comparison with the experimental data, we have to use some reasonable, simple and 
sufficiently flexible parameterization for the unknown functions //y and Vq. We suggest 
the following. 

a) Normalized distribution fiy is assumed in the form 



Pv(%o,Q 2 ) 



r(r + s + 2) 

r(r + i)r(s + r 



x r (i - x y 



< Xq < 1, 



(3.32) 



where the Q 2 — dependence is involved in the parameters r, s. 

b) For the function Vq{x)x we shall use the parameterization suggested in Eqs. (Q2.34)- 
(Q2.38) 



Vq(x)x — ^fCnorm • f(x, Xq), f(x, Xq) = 



:i - x) i - 



X 



a/2x 



(3.33) 



where the upper (lower) sign in the first relation refers to the region x > x (x < x ) 
and 



f{x,x ) 1 



x 'o 

X 2 



n -1 



dx 



(3.34) 



which follows e.g. from Eqs. (Q2.34), (P3.25). Now, apparently one has to accept the 
parameter a ~ raj (E kin ) depends on x as well. Let us consider the following. Sequence 
of the pictures in Fig. |l]a can be understood as the pictures of the one and the same 
nucleon, but "seen with different resolutions" as outlined in Fig. [l]&. Then, it is natural 
to assume the four-momentum P of the parton from some picture can be obtained from 
four- momenta p\ of n partons in a picture more rightwards, representing the parton "seen 
with better resolution" : 



(3.35) 
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Obviously, the mean values satisfy 

(P ) = n (pxo) , ( P > = Ccopp ■ n (|pa|) , < c corr < 1, (3.36) 



where c corr equals 0(1) for the extreme case, when the motion of the partons in the 
corresponding subset is completely uncorrelated (correlated). The last relations imply 
the effective masses and kinetic energies obey 

m{P) > n ■ m(p), (E kin (P)) < n ■ (E kin (p)) , (3.37) 

which means the quantity a is a non- decreasing function of Xq. In this moment we know 
nothing more about this function, in the next section we shall show a reasonable agreement 
with the experimental data can be obtained with 

a (x ) = ai (- ln(xo)r° 2 • (3-38) 

Since we parameterize the function V Q ' rather than the function p, it will be useful the 
quantity e(xo), defined in Eq. Q3.8|) and afterwards entering the important Eq. ( p. 171) , to 
express also in terms of Vq. Obviously, using Eqs. ( |3.11| ) and ( |3.33| ) one gets 



e(x 



P(e, x )ede = - £ V&)t (t + | ) f (l - 

Mr 1 ( x 4 \ 

= c norm — J /(£, x ) - f) d£. (3.39) 

Now, we can our results shortly summarize. If there are given some values of the 
free parameters c q , a, r, s, a±, at2, then the corresponding proton and neutron structure 
functions can be directly calculated according to Eqs. ( |3.16|) , ( 3.171) , ( |3.31|) , where the 



distribution \xy is given by Eq. (|3.32| ), the function e(xo) by Eq. (|3.39|) with the use of 



Eqs. ( p.33| ), (|3.34| ), ( |3.38|) and the functions (p,tpi,ip2 are: 

<p(x, xq) = 2x 2 [ ri(Z,x )d£, (3.40) 



ifi 1 (x,x )=x I r](£,x ) x l — % l -d£, (3.41) 

Mx,x ) = f f\{^)-F== ~ 2 2< (3.42) 



where 



Xq) = Cnarm0{£ ~ Xq) ^ ' ° 2 ■ (3.43) 



e+xi 

The last expression is calculated from Eqs. ( p.33j ), ( |3.34| ) and ( |3.38| ) with the use of the 
step function 9(y) = +1(— 1) for y > (y < 0). 
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4 Comparison with the experimental data 



Now we shall try to compare our formulas for the structure functions with the existing 
data. We shall not attempt to make a consistent, global fit of the free parameters based on 
some rigorous fitting procedure, but only show the set of optimal parameters obtained by 
their tentative varying on the computer "by hand". Moreover, our constraint will be only 
agreement with the proton structure functions F 2 and g\. It means that the parameter a, 
controlling asymmetry between the proton and neutron spin functions, must be somehow 
fixed before the fitting. For the first approximation we use the SU(6) value, a = 2/3 [see 
item 6c) in the previous section]. 

For a comparison with F 2 we use the parameterizations of the world data suggested 
15] and |16|], both taken for Q 2 = lOGeV 2 . The data for g\ are taken over from the 



in 



recent paper [16 of the SMC Collaboration. After some checking on the computer, the 



optimal set of the free parameters is considered: 

c q = 0.43, r = -0.49, s = 6.5, ot x = 1.6, a 2 = 1.5 (4.1) 

Results of the calculation of the proton structure functions g\ and F 2 with these parame- 
ters are shown in Figs. [| [| together with the data. Let us remark, the experimental points 
for gi correspond to the values evolved in [|TB| to Q 2 = lOGeV 2 . In the error bars all the 



quoted errors (statistical, systematic and those due to uncertainty of QCD evolution) are 
combined. Obviously, the agreement with the experimental data in both the figures can 
be considered very good, particularly if we take into account that our parameterization of 
the unknown distributions is perhaps the simplest possible and moreover, the parameters 
( |4. 1| ) still may not be optimal. 



Now, having "tuned" the free parameters by the g\ and F 2 , one can predict the re- 
maining functions g 2 and F 2va i. The results are shown in Figs. |4], [5]. Our xg 2 surely 



does not contradict the experimental data [17], which are compatible with zero - with 



statistical errors substantially bigger, than the vertical range of the figure. But instead 
of the data, the comparison is done with Wandzura Wilczek |18|] twist-2 term for xg 2 vw , 
which is evaluated in |17| from the corresponding g\. It is obvious, that two completely 



different approaches give at least qualitatively very similar results. The proton valence 
function F 2va i in Fig. [5] is compared with the corresponding combination of the distri- 
butions xuy(x) and xdy(x) obtained (for Q 2 = AGeV) in |]TD| by the standard global 
analysis: 

F 2va i(x) = -xuy(x) + -xdy(x), J uy(x)dx = J dy(x)dx = 1. (4.2) 

Apparently, the agreement can be considered good. Nevertheless, let us note, a part of 
the difference between both the curves can be due to our simplification: uy ~ dy. On the 
other hand, the distributions of the u- and (i-valence quarks are not equal and enter the 
functions g\ and F 2vai with different weights, see Eqs. ( |3.22| ), (fO|). At the same time, 
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our prediction for F 2va i is based on the combination valid for g± 



a^u v (x) + (1 - 2a)^d v (x) w 1 + 6a qv(x), (4.3) 

which is related to F 2va i 

8 1 9/8 1 

F 2 vai(x) w + -xq v (x) = xq v « - - — I -iuv(i) + -sdy(a;) 

16 1 

= — xu v (x) - —xd v (x) (4.4) 
15 15 

for our value a = 2/3. The last combination is also shown in the Fig. [| and apparently 
gives a slightly better agreement with our simplified F 2va i- In any case, one can note, that 
the two different procedures, the standard one (uses input on F 2 , Fg + QCD) and our 
(uses input on F 2 ,gi + our model) give a very similar picture of the function F 2va i{x), 
which is not directly measurable. 

5 Discussion 

Let us make a few comments on the obtained results. First of all, it should be pointed 
out, that our structure functions in Figs. are calculated on the basis of very simple 
parameterization of the unknown distributions fi(xo) and Vq(x, Xq), but on the other 
hand it is essential, that the contributions from the individual components of the quark 
distribution correctly take into account the intrinsic quark motion, which is particularly 
important for the spin structure function. The effect of this motion on g\ is demonstrated 
in Fig. Ql and the fact, that we succeeded to achieve a good agreement with the data also 
in Fig. is just thanks to this effect. For a better insight, how our structure functions 
are generated, in the Fig. || we have displayed the initial distribution function Vq(x, xq) 
drawn for a few values xo, together with the corresponding structure functions F 2 ,gi, xg 2 . 
The complete structure functions are their superpositions - weighted in the corresponding 
way with the use of the distribution fiv(xo)- 

Further, also some other assumptions of the model are possibly oversimplified, for 
a more precise calculation, at least some of them could be rightly modified - but at 
a price of introducing the additional free parameters. For example, the constant w sea 
[see Eq. (|3.13|)] should take into account some suppression of the s— quarks [O] and 



probably should allow a weak dependence on x . Also for the constant c q [see Eq. (|3.5|)1 
some Xo— dependence should be allowed. Concerning this constant, let us make one more 
comment. The standard global fit |19| suggests (at Q 2 = lOGeV 2 ) the quarks carry ~ 56% 



of the nucleon energy and our fitted value c q from the Eq. (|3.5|) is roughly 43%. This 
difference is mainly due to the different relations between the distribution and structure 
function in both the approaches, see Eqs. (P3.38), (P3.59). The second relation (valid for 
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a subset of quarks with effective mass Xq), multiplied by x 2 and then integrated by parts 
gives 

/ xF(x,x )dx = \ j F 2 (x, x Q ) { 3 + ^ ) dx, (5.1) 
Jx% 4J X 2 \ x 2 J 

which for the static quarks [F(x, xq) ~ F 2 (x, xq)/x — 5(x — xo), see discussion after Eq. 
(P3.59)] coincides with the standard relation. Nevertheless, generally both the relations 
imply different rate of the nucleon energy carried by quarks. One can check numerically 
that for our F 2 (x, x$) in a dominant region of x$ the term (x/xq) 2 in the integral ( |5.1|) plays 
a minor role (see also Fig. |6], positions of the maxima of F 2 's are above the corresponding 
Xo, in particular for lower x )> so as a result we get 3/4 of the standard estimation of the 
quark contribution to the nucleon energy. This ratio agrees with the ratio obtained from 
the corresponding fits: 3/4 ~ 43% / 56%. 

In the previous section we mentioned the effect of different shape of the u— and 
d— quark distributions. Also the result g™(x) = following from Eq. ( |3.23| ) is a conse- 
quence of the simplifying assumptions u n = d n = u p = d p = h and a = 2/3. A more 
general approach u ^ d would give 

H n ( P o) = ^ (-u n ( Po ) + d n (p )) = 1 (-d^ipo) + u p (p )) , (5.2) 

so u 7^ d implies H n ^ 0. In fact, the global fit analysis proves, that d p val (x) is slightly 
"narrower" than u p val {x). It means, considering very qualitatively, in accordance with the 
last equation, in the function g™{x) the negative term should dominate for smaller x, 
which does not contradict the experimental data. Obviously, a proper accounting for 
the difference u ^ d into the model should enable to calculate consistently in a better 
approximation not only the proton and neutron structure functions F 2 ,gi,g 2 , but also 
the neutrino structure functions. Apparently then one could make a "super-global" fit 
covering the both unpolarized and polarized DIS data. As a result, the flavor- dependent 
quark distributions V (x,x ) [or equivalently p(e,x )] together with the corresponding 
effective mass distributions and the parameter a controlling the relative spin contribution 
of the u— and d— quarks, could be obtained. 

Finally, let us point out, inclusion the spin structure function into the fit in our model 
enables to obtain some information about the distribution of the quark effective masses. 
Within our approach there are two distributions, /iy and fi, relevant for the description of 
the quark effective masses in the nucleon. The extrapolation of our parameterization for 
the \i distribution with the use of the relations a ~ raj (E kin ) and ( |3.9| ),( p.32| ), (|3.38|) ,(|4~T 
give for x ~ ¥ : 



, . M^o) x X Q , ~ ,, , 

M*o) ~ J7-T -> 77Z ,, F v -> 7^ 7U, ( 5 - 3 ) 



which implies the extrapolated \x is not integrable in this limit. On the other hand, the 
basic distribution /iy, parameterized by Eq. ( p.32|) with the r, s from the set ( |4.1|) and 
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with the use of the known relation zT(z) = T(z + 1) can give an estimate of the mean 
value: 

(x ) v = / + + a l 2 * 0.064, (5.4) 

i.e. (m) ~ 60MeV for Q 2 = lOGeV 2 . The corresponding kinetic term calculated as 

(E k in) v = J fi v (x )(e(x ) - Mx )dx (5.5) 

gives a similar number (~ 60MeV). Let us recall, these numbers are related to the valence 
quark term. The Q 2 -dependence is involved only in the distribution fi v (xo,Q 2 ), i.e. in 
our parameterization ( |3.32| ) only via the powers r(Q 2 ), s(Q 2 ). It follows, the structure 
functions, which enhance in a low— x region for increasing Q 2 , must be generated by the 
distribution nv(xo, Q 2 ) in which the mean effective mass {xq) v drops for increasing Q 2 - 
in the qualitative agreement with an intuitive expectation. 



6 Summary 

In the present paper, with the use of the results obtained in the preceding ones we 
proposed an alternative covariant formulation of the QPM. The initial postulates of the 
standard and our approach are basically the same, despite of that the relations between 
the structure and distribution functions obtained in both the approaches are not identical. 
It is due to the fact, that in the standard approach the intrinsic quark motion is effectively 
suppressed by the use of the approximation = xP^. On the other hand, we have shown 
the master equations can be solved without the use of the this approximation, so in the 
corresponding solution the quark intrinsic motion is consistently taken into account. On 
the basis of the obtained relations (a priori valid for the version of naive QPM - with non- 
static quarks on mass shell) we propose the model, in which the distributions (/j,v, V ) 
reflecting the parton dynamics are introduced with some free parameters. With the use 
of this model we calculated simultaneously the proton structure functions F 2 , F2vai > Qi > 92, 
assuming only the valence quarks term contributes to the proton spin. By a comparison 
with the data on F 2 ,gi {Q 2 — WGeV 2 ) we fix the free parameters and then on the basis 
of these parameters the functions F 2va i and g 2 are predicted. We found out: 

1) Both the unpolarized structure functions are well reproduced by the model. The 
comparison is done with the data on F 2 and with the F 2va i obtained from the standard 
global analysis data. 

2) At the same time, the model well agrees with the data on g\. The calculated g 2 
does not contradict the existing experimental data as well, but since in the data there are 
still rather big statistical errors, it is hard to say more. 

3) Analysis of the fixed parameters within our approach suggests: 

i) The quarks carry less the proton energy (almost by the factor 3/4), than estimated 
from the standard analysis. 
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ii) The average effective mass related to the valence quark term can be roughly 60MeV 
and a similar energy can be ascribed to the corresponding motion. 

So on the end we can underline the proposed model offers apart of the other results 
also the consistent explanation, why the experimentally measured proton spin function gi 
is less, than it is predicted in the standard approach based on the QCD improved QPM. 



7 Appendix 

In this appendix we shall explain in more detail, that approximation (|1.1|) is sufficiently 
correct for the purpose of this paper. Structure functions given by Eqs. (P3.52), (P3.47) 
and (Q2.13), (Q2.14) are given by the integrals of the form 

/( PQ \ Q 2 
K(Po,Pi)S [j^ ~ x ) d 3 P, Pl=pi + pl+pl + m 2 , x= 2Mv' ^'^ 

In accordance with Eqs. (P2.8)-(P2.13) one can express 



pq Po+Pi^Pi Pt //4M 2 M 2 \ 2M 

+ T X -T7\ / ~™ IT ) x —xcosip (7.2) 



I(x,a) = J K{p ,p 1 )8 ~ a Pr cos (p - x^j dp x p T dp T dip, (7.4) 



Mv M k My V Q 2 K J 

and if one assumes ko 3> M, then the second term can be neglected and after introducing 
the variables pt, 

P2 = Pt cos (p, p 3 = p T sin <p; d 3 p = dpip T dp T dtp (7.3) 
the integral (|7.1| ) reads 



where 



M V V Q 2 k o J k > 

Obviously, our approximation ( |1.1|) used in Eqs. (P3.53) and (Q2.15) corresponds to 
a = i.e. Q 2 — ► oo and now we shall try to estimate the effect of the finite Q 2 in Eq. 
( [T.4|) . First, we shall integrate over p\ using the rule 

5(f(x))dx = J2 H f,7 X ^ dx, f(x j )=0, (7.6) 
j J \ x jl 

which with (just one p{) root in the argument of the 5— function in Eq. ( [T.4| ) gives 

P0 + Pl \ 7 PO X f V , /- -S 

— apT cos(p — x dpi = o {pi — pi) dpi, (7.7) 



r^-^ ™. (7.5) 



M J x + a PT cos ip 
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where 



Po = - 



M (x + apr cos <p) + 



Pt 



M (x + apr cos 99) 



(7- 



1 

Pl = 2 



, , , . m 2 + p% 

M [x + apr cos ip) 



(7.9) 



M (x + apr cos (p) 
Now the integral ( |7.4| ) can be expressed 

J(z,o)= / X(p ,Pi) — ■ — — p T dp T dp. (7.10) 

J x + apr cos (p 

The upper limit of can be with the use of Eq. (P3.ll) calculated from the equation 

PO = ^max = — (1 + X 2 Q ) , X = Ul/M, (7.11) 

which implies the equation 

Ap 2 T + Bp T + C = 0, (7.12) 

where 



A = 1 + MVcosV B = 2Macosp(xM - E max ), C = -M 2 (l - x)(x - x 2 



(7.13) 



and correspondingly 



PTmax = ^4 • ( 7 - M ) 

Since one requires px > and obviously we have A > 0, C < in our region x\ < x < 1, 
so only the upper sign root is acceptable. For Q 2 — > oo this root reads 



PTmax = Myj{\-X){x-Xl) (7.15) 

in accordance with Eq. (P2.19). The integral (|7.10|) can be now expressed 



(*1tv rpT max (x,acosip) 

I(x,a) — I / L (x + apx cos (p,pt) dpxdcp, (7-16) 



o jo 



where 



L (x + ap T cos (p, p T ) = if(p ,Pi) — , ° T • (7.17) 

x + ap T cos 
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Now, using the rule for differentiation of the integral 



ffl = J'® h(z,t)dz^f t = J 



9{t) dh(z,t) , , , ,dg 

-^dz + h(g(t),t)f t (7.18) 



and taking into account our integrated function vanishes at the boundary of phase space, 
L (x + ap Tmax cos <p,Pt max) = 0, then the second term in the last equation equals zero and 
the integral ( |7.l6|) can be expanded 



I(x, a) 



E 

j=0 



d j I(x, a) 



da? 



(7.19) 



a=0 



where 



d j I(x, a) 



da? 



a=0 



2w rvT^Aco^) giLjx + ap T cos y,p T ) 
o Jo 



da? 



dp T d(p. 



(7.20) 



a=0 



At the same time, one can write 



(x,acosip) Qj 

Jo 



dprdip 



a=0 



2tt np T m ax (a, a COS tp) Qj 

o Jo 



dxi 



(pt cos ifY dprdip 



a=0 



(7.21) 



i.e. the j— odd terms vanish, so the difference I(x,a) — l(x,0) can be estimated by an 
expression of the magnitude of the second order in the parameter a. Quantitatively, this 
magnitude depends on the actual shape of the function K. So let us now calculate the 
integral ( |7.l6| ) for the functions Vo,F 2 ,gi,g2 in accordance with the parameterization 
used in this and the previous paper. First of all, one has to replace the parameterized 
distribution function Vo(x)dx by the corresponding G(po)d 3 p. The equations (Q2.34), 
(Q2.37) and (Q2.38) imply 



P(Po) 



which gives 
P(po) ~ 



I, 2 2 Po 



X) 



a/2x Q 



1 _ Z£ 

x 



a/ 2xo 



Po 



Po 



a/2x 



M I Xr, 

— (l + — 

X 



M , 2 , 

-C'max — (1 T Xq J . 



(7.22) 



Normalization factor is calculated by integration 



N 



Po 



a/2x 



dpo 



2 a/2x + 1 V l+^o 



M 1 +xl 



2 \ a/2xo+l 



(7.23) 



(7.24) 
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The distribution P(p )dpo is related to the distribution G{po)d 3 p by Eq. (P3.14), so we 
obtain 



Anpsm^J (pq/ m) 2 — 1 

where a is for the numerical calculation parameterized in accordance with Eqs. (|3.38| ) and 
( fill) . Now, the integrals corresponding to the functions Vo, F 2 , gi, g 2 could be calculated 
without approximation ( |1.1| ). In an accordance with Eqs. (P3.47), (P3.52), (Q2.17), 
(Q2.18) one has to compute 

V (x, x ) = [ G(p )5 (£L - x) d 3 p, (7.26) 



Mv 

F 2 (x,x ) = x 2 M J G(p )6 (jfc-x) ^f, (7-27) 

9l (x,x )=w spm - J G( Po )-j===5(—-x) — , (7.28) 

7TI f Pi / PQ \ d^p 

g 2 (x,x ) = -w spm - J G(p ) \m^- x ) — > ( 7 - 29 ) 



where w sp i n is the factor given by Eq. ( |3.24| ). The last integrals are the special cases of 



the integral (|7.1|) , so they can be rewritten in the form of Eq. ( [7.1 6|) allowing an easy 
numerical calculation 

V {x,x )= / g(po) ■ P ° PT dprdtp, (7.30) 

Jo Jo x + a PT COS if 

p2tt rPTxaax{x,acOBip) 

F 2 (x,x )=x 2 M / G(p ) , dprdip, (7.31) 

Jo io ^ + a PT COS y? 

#i(x,x ) = w spin — / / G(p )— === — ■ dp T dLp, (7.32) 

2 Jo Jo VPo - Pi x + a ^ T cos V 

g 2 (x,x ) = -w spin — / / G(p )^=== — ■ dp T dip, (7.33) 

2 Jo Jo y/p 2 -p 2 x + ap T COS (f 



20 



where quantities G, a, p , pi, PTmax are given by Eqs. ( |7.25| ), (\H$, (|7^), <JT% , ( |7.14| ) 



respectively. These functions F 2 ,g\,g2 could replace the functions (|3.11|) , (|3.24| ) used in 



the present paper. In the Fig. [7| we have shown the results of numerical calculation of the 
functions ( 7.30 ) - (|7.33| ) for a = 0, which is equivalent to the approximation (LI) and for 



a > 0, corresponding to Q 2 = 3GeV 2 /c 2 and k — > oo. The curves a = also correspond 
to the results of calculation in Fig. phased on Eqs. ( |3.33j ) and ( |3.40|) - ( |3.43|) . The figures 



suggest convincingly, that the effect of finite Q 2 even for its rather low value, can play 
some role in our integrals ( |7.30| ) - ( [7.33j ) only for the higher effective masses (xq ~ 0.5), 
which are in our approach evidently on the tail of the mass distribution /iy(x )- It means, 
that the approximation ( |1 . 1| ) is for the purpose of our model fully acceptable. 
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Figure 1: Nucleoli consisting of the valence and sea quarks - with different resolutions, 
see text. 




Figure 2: Proton spin structure function gi(x) at Q 2 = lOGeV 2 . The points represent 



experimental data the curve is the result of our calculation. 
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Figure 3: Proton structure function F 2 (x) at Q 2 = lOGeV 2 . The dotted and dashed 
curves represent the fits of the experimental data suggested in j0| and [|T^]. The full 
curve is the result of our calculation. 




Figure 4: Proton structure function xg 2 (x) at Q 2 = lOGeV 2 . The points represent the 
term xg 2 vw (x) (see text), the curve is the result of our calculation. 
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Figure 5: Proton structure function F2 Va i(x). The dashed and dotted curves represent the 
functions based on the standard global analysis according to the relations ( |4.2| ) and ( O ) , 
see text. The full curve is the result of our calculation. 
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Figure 6: Distribution functions V (x, x ) drawn for x = 0.005,0.015,0.05,0.15,0.5 and 
the structure functions generated correspondingly. The calculation is based on the Eqs. 
(PI), ( ETTTD , (IPD and flOOp . 
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Figure 7: Distribution functions Vq(x, Xq) and the corresponding structure functions 
drawn for x = 0.005,0.015,0.05,0.15,0.5. The functions are calculated according to 
Eqs. (pD - (|7^3| ) for Q 2 -> oo (full lines) and Q 2 = 3GeV 2 /c 2 (dotted lines). 
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